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In this paper, a study on modeling of a piezoelectric transducer type rotary 
traveling wave ultrasonic motor (USM) is presented. First a mathematical 
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model and numerical simulation results are achieved. The model is based on 
the theory of piezoelectricity and physic theory. An experimental model is 
worked out and compared to the numerical model. The influence of the 


temperature on characteristics such as the rotational speed of the motor is 

considered. The speed of the USM is measured at temperature between 17°C 
Keywords: and 50°C. To develop suitable control strategies for the drive, a fuzzy model 
type Takagi-Sugeno is used. The unknown parameters of the output 
membership functions are determined by least square method. Experimental 
data are used to examine the validity of the fuzzy model. Comparison 
Temperature effect between experimental and calculated data of the fuzzy model indicates that 
Ultrasonic motor the fuzzy model can well describe the nonlinear characteristics among the 
USR60 frequency of driving voltage and rotating speed. 
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1. INTRODUCTION 

Ultrasonic motors are a type of actuator which have some excellent performances such as high 
holding torque, high torque at low speed, quiet operation, simple structure, compact size and no 
electromagnetic interference [1]. 

Due to such characteristics, ultrasonic motors have been in many practical applications, such as in 
robots, medical instruments, cameras, aeronautics, MEMS and many others. Precise theoretical modeling of 
motion and energy conversion in ultrasonic motors is of great interest both from theoretical and practical 
viewpoints. One of the most interesting types of motors is a travelling wave ultrasonic motor. However, it is 
difficult to derive complex mathematical model of the USM. Moreover, the control characteristics of the 
motor are highly nonlinear. 

The exact values of motor parameters are difficult to obtain and they are time-varying due to 
increase in temperature and changes in motor drive operating conditions such as driving frequency, and load 
torque. Several attempts have been carried out to develop the model of the motor. The equivalent circuit was 
widely used on ultrasonic motor modeling [1]-[6]. But this type of model is a rude modeling method of the 
dynamic characteristics of ultrasonic motor and cannot be used to design motor’s controller directly. 
Dongsheng Zhang, Shiyu Wang and Jie Xiud [7] developed an analytical dynamic model of ring ultrasonic 
motor. Based on this model, the piezoelectric parametric effects on the wave distortion and contact 
mechanics are examined. In [8]-[10] differential evolution algorithm is used for motor’s model identification. 
The algorithm reduces the time required in the process of identification. Many researchers considered 
mechanical characteristic of contact frictional layer and calculated displacement and speed of the rotor using 
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finite element method [11]-[13]. In [14]-[16] dynamic contact problems between a rigid rotor and deformable 
stator were approximately solved. Due to the limitations of theoretical modeling, control design usually 
adopts identification method based on the tested data. According to the different identification methods, the 
forms of model can be the transfer function, differential equation, neural network, and so forth. In recent 
years, fuzzy modeling method based on fuzzy reasoning is gradually arisen. The same as in the neural 
network model, fuzzy model is also based on experimental data, easy to show the nonlinear information. The 
fuzzy method is mostly used to realize speed and position control [17]-[19]. 

This paper presents both a numerical model and a fuzzy modeling method based on Takagi-Sugeno 
reasoning is used to obtain the model of ultrasonic motor. The Takagi-Sugeno model is chosen to make the 
proposed model easy to use on control design. The output of Takagi-Sugeno model is a linear mathematical 
expression. Contrariwise, the Mamdani model gives a fuzzy output. The proposed work considers the effect 
of temperature. This work will be very useful for the optimal design, reducing the heat loss, preserving the 
desired rotary speed, and improvement of precise control of USM. This method is based on fuzzy model 
structure and experimental measurements. The data contains the driving frequency, the temperature and 
rotary speed of ultrasonic motor. The least squares method is used to identify the unknown parameters of the 
conclusions of fuzzy rules. 

This paper is organized as in Section 2, a mathematical model is presented. Simulations are 
achieved on MATLAB. In Section 3, experiments are performed to verify the effectiveness of the analytical 
model. The effect of temperature is discussed. An adaptive fuzzy model is investigated in Section 4, taking 
into account the motor’s non-linearity. Finally, conclusions of this work are presented in Section 5. 


Table 1. Design specifications and simulation parameters of USR60 


Items Descriptions Values 

f Supply frequency 40-42kHz 

u Supply voltage 100V ims 
Speed Rated speed 90rpm 
Power Rated power 3w 

M, Load torque 0-0.32Nm 
M Stator mass 10.1e-3Kg 
MR Rotor mass 30e-3Kg 
dz, dp Rotor damping 5% 

D Stator damping 15.4Ns/m 

C Stiffness 5.95e+8N/m 
Fy Applied force 160N 

J Moment of inertia of the rotor 7.2e-6K gm? 
0 Coupling factor 0.2263N/V 


2. USM MODELING 

In [20] an equivalent mechanical two-mode approximation is given, representing the modal 
amplitudes w,, and wz; of the sine- and cosine-mode of the vibrating system (stator and ceramic), see 
Figure 1. The parameter m describes the modal mass of stator and ceramic, and dy, represents the structural 
damping. Cs and cç reflect the equivalent mechanical stiffness of the stator and ceramic, respectively. y4 and 
Y, are displacements which are proportional to the electric currents flowing into the ceramic layer. Cross- 
coupling between the two modes is represented by the small disturbances €,, €2. The determination of the 


model parameters can be based on energy considerations or detailed finite element analysis of the stator [21]- 
[24]. 








Figure 1. Mechanical model for the (stator + ceramic) of USM 
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2.1. Mathematical model 

An ultrasonic motor uses mechanical vibrations in the ultrasonic range as its drive source. A voltage 
is applied to the piezoelectric ceramic elements through two electrodes to generate alternating expansions and 
contractions in the ceramic body. The magnitude of these oscillations is of the order of lum. In order to 
obtain a higher gain; the resonance effect of the ceramics in the ultrasonic range is used. The specifications of 
the USR60-object of this study-are given in Table 1. The stator of USM is considered as a piezoelectric plate, 
polarized in its thickness direction [1]. 

When two specific phase voltages are applied to the piezoelectric elements, a travelling wave is 
generated on the surface of the stator. The travelling wave produces vibration with an elliptic locus on the 
surface of the stator. This stator vibration is transformed into a rotation through friction contact between the 
stator and rotor. The rotor is pressed against the stator by a disc spring (160N). The contact area between 
stator and rotor moves with the travelling wave. The rotor is modeled by a rigid strip [1]. In dealing with the 
dynamics of the rotor, two degrees of freedom must be taken into account: first the rotation of the rotor and 
second the motion in Z direction, as shown by [25]. The dynamics of the vertical rotor motion is obtained 
through the force equilibrium in Z direction. 


MpZ+dzz = F; — Fy (1) 


with mp is the mass of the rotor, and dz is an appropriate damping term of the vertical motion. The applied 
axial load, Fy, is a system input, and Fz is the force at the contact area. The equation of rotational motion is 
calculated by. 


J6+dp0 = Mp — M, (2) 


where M, is the applied torque (mechanical load), dp denotes the damping in spinning direction and J 
represents the moment of inertia of the rotor. The Mp term is the frictional pressure at the contact area. The 
dp term is small, but not negligible. The dynamics of the stator can be described by 


M-:w,+D-w,z+C:w,=O0-u (3) 


where the parameter M describes the model mass vector of stator and ceramic, D reflects the structural 
damping vector, and C represents the equivalent mechanical stiffness vector of the stator and ceramic. The 
vector u is the phase voltages, and O represents the electromechanical coupling matrix considering the cross 
excitation. In order to consider the interaction between the stator/rotor-contact, Equation (3) reflecting the 
stator’s dynamics is to extend by a nonlinear modal force vector F4 yielding. 


M-:w~.+D-w,z+C:w, =O-+ut Fy (4) 


Hence, the dynamics of the motor is completed. An integration of this modeling approach based on closed 
form solutions for the contact mechanism into the simulation model is feasible. 


2.2. Simulation results 

The implementation is performed on MATLAB and the parameters of USM are scheduled in 
Table 1. The temperature of the motor is supposed to be 20°C. The speed of the motor has its maximum at 
the mechanical resonant frequency. It is due to the fact that the revolving speed of the motor is proportional 
to the vibration force of piezoelectric elements. So, any deviation from this frequency degrades the motor 
performance. However, this effect seems more serious for frequency decrements. 

To avoid these phenomena, speed of USM was controlled in 40-42kHz frequency range. Figure 2 
shows speed-frequency characteristic of USM under different load torques. Speed-frequency characteristic of 
USM is not linear. Rotary speed decreases as the load torque increases. As mentioned, variations in driving 
frequency changes speed of USM. Figure 3 represents speed-torque characteristic of USM for driving 
frequency f = 40kHz. When the load torque exceeds 0.32Nm, the model becomes imprecise. 
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Figure 2. Speed-frequency characteristics of USM Figure 3. Speed-torque characteristic USM 


3. EXPERIMENTAL DATA AND RESULTS 

To obtain the temperature and the performance characteristics of the motor, an experiment setup is 
established. The motor is alimented by the two-phase voltages with 90° phase shift. The teensy development 
board running at 72MHz and the driver circuit are utilized to generate, amplifier, and filter the applied 
voltages of the USM. 

A temperature sensor (LM35) is used to measure the temperature of the stator in the temperature 
range of 17°C to 50°C. A highest speed of the USM can be obtained when the stator is excited by the 
resonant frequency. The driving frequency is measured by oscilloscope. A capacitive encoder (AMT102 with 
1024 pulses per revolution) is used to measure the rotating speed of USM. The development board is used to 
measure data from both the position and temperature sensors. The experimental bench and its block diagram 
are shown in Figure 4. 

A comparison between simulation and experimental results for speed-frequency characteristic of 
USM under different load torques, is given in Figure 5. The figure shows that the simulation results are 
compatible with the experimental values. However, for frequencies above 41.1kHz we notice a slight 
deviation due to the increase in motor temperature. The matlab model is not time dependent. Nevertheless, 
during motor operation and as time passes, the temperature of the ceramic elements increases, which 
decreases the USM speed. 

In order to measure the temperature of the USM, the stator is excited by the resonant frequency. 
Figure 6 shows the speed of USM versus time. It is evident that speed decreases according to the time 
dependent temperature rises. Figure 7 shows the measured temperature versus time. After 13 minutes, the 
motor temperature exceeds 50°C. The temperature-time curve is almost linear. In terms of its effect on 
resonant frequency, however, temperature change is a more important factor than dimensional precision. The 
effect of temperature is shown in Figure 8. From this figure, it is concluded that the resonant frequency- 
temperature characteristic is linear from 20°C to 50°C, and it can be approximated by a linear equation [26]- 
[27]. In [1] the authors affirmed that this characteristic is linear from -10°C to 55°C. 
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Figure 4. Experimental bench and its block diagram 
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Figure 5. Comparison between simulation and experiment results 
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Figure 8. Temperature-dependence of resonant frequency 


4. ADAPTIVE FUZZY MODEL 

The ultrasonic motor has nonlinear speed characteristics which vary with drive operating conditions 
such as temperature rise and applied load torque changes. Furthermore, the motor has no mathematic model 
for control by far because the drive principle is different from the general electromagnetic type motors. In 
this section, a fuzzy model is proposed to solve these problems from the point of utility. The input of the 
fuzzy model is the frequency f, and the output is the speed w. The frequency f is defined in the range from 
40 to 42kHz. Looking at the input-output data shown in Figure 5, we can see that the speed-frequency 
characteristics change as the frequency decreases. Five fuzzy sets: “very_small”, “small”, “medium”, “big”, 
and “very_big” are defined as shown in Figure 9. We have a model with the five implications. 


if f is very small then œw, = a'f + bi (5) 
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if f is small then w, = a'f +b, (6) 
if f is medium then w3 = a; : f + b} (7) 
if f isbig then wœ, = a,'f + b, (8) 
if f is very_big then ws = as 'f + bs; (9) 


where w; (1= 1...5) is the rotary speed. The five membership functions of Figure 9 are given by. 


m, (f) = exp(—o(f — 40000)?) (10) 
m,(f) = exp(—o(f — 40500)?) (11) 
m3(f) = exp(—o(f — 41000)?) (12) 
m,(f) = exp(—o(f — 41500)°) (13) 
ms(f) = exp(—o(f — 42000)?) (14) 


with ø = 6,65e~°. Then the motor speed is obtained by defuzzification. 
w = Xi pi’ @ = Pry (15) 


with B; =m,;/X2m,;, PT = [a,b], and Y = [f,1]". In this paper the vector P is calculated by the least 
squares algorithm. 
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Figure 9. Membership functions of fuzzy sets 


Figure 10 shows the experimental model (no load) and the fuzzy model of the speed-frequency 
characteristic. Such characteristic indicates the similarities between the experimental nonlinear model and the 
fuzzy model. This model can then be used to perform a direct control design. With the adaptive fuzzy 
approach, motor’s nonlinearities are masked by local linearization using the Takagi-Sugeno reasoning. The 
same conclusions are valuable for load below 0.3Nm. 
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Figure 10. The speed-frequency characteristic 


5. CONCLUSION 

In this paper, a study on modeling of USM is presented. A numerical model of the motor is 
achieved. The simulation results are compared with the experimental ones which give good agreement of the 
model. The effect of the temperature on the characteristics of the USM is discussed. Furthermore, the effects 
of driving frequency and load torque are identified. Based on the experimental data, a fuzzy model is 
achieved. The model can well simulate the nonlinear relationship among the driving frequency and rotary 
speed. Results of the experiments and measurements showed the validation of the fuzzy modeling. The paper 
shows that the nonlinear characteristics of USM are fit to be described. Furthermore, the presented method is 
general, and can be extended to other travelling wave type motors. To make precise control of USM, the 
proposed fuzzy model can be used. A fuzzy PID regulator is feasible to balance the variations in the 
Operating conditions and preserve the desired rotary speed. As a future work, fuzzy logic will be combined to 
sliding mode technique to achieve a robust control scheme. The main advantage of fuzzy sliding mode 
control is that the control method achieves asymptotic stability of the closed-loop system. 
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